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We describe an efficient practical procedure for enumerating and regrouping vacuum Feynman 
graphs of a given order in perturbation theory. The method is based on a combination 
of Schwinger-Dyson equations and the two-particle-irreducible ("skeleton") expansion. The 
regrouping leads to skeletons containing only free propagators, together with "ring diagrams" 
containing all the self-energy insertions. As a consequence, relatively few diagrams need to 
be drawn and integrations carried out at any single stage of the computation and, in low 
dimensions, overlapping ultraviolet /infrared sub divergences can be cleanly isolated. As an 
illustration we enumerate the graphs contributing to the 4-loop free energy in QCD, explicitly 
in a continuum and more compactly in a lattice regularization. 
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1. Introduction 



There are many physics contexts where multiloop Feynman diagram computations are carried 
out. In QED one goes up to 4-loop level (for reviews see, e.g., Q), because experiments are 
so precise. In particle physics phenomenology, particularly QCD, one goes up to 4-loop level 
(see, e.g., Q), because the coupling constant is not small. In studying critical phenomena 
in the simplest 0(N) condensed matter systems, one goes up to 5-loop level (see, e.g., fH), 
because the effective expansion parameter is not small. 

Studies of QCD at a finite temperature T are faced with a similar challenge. Indeed, the 
coupling constant expansion converges even worse than at zero temperature, requiring at least 
T ^> 10 3 Aqcd to make any sense at all [|], HJ. So far, though, only resummed 3-loop level has 
been reached for the simplest physical observable, the free energy ||, because a broken Lorentz 
symmetry makes the analysis much more complicated than in the cases mentioned above. In 
fact, even in principle only one more order is (partly) computable, and then the expansion 
breaks down completely |7|]. Multiloop computations are not useless, though: the infrared 
problems can be isolated to a simple three-dimensional (3d) effective field theory || and 
studied non-perturbatively there Q, but to convert the results to physical units from lattice 
regularization still necessitates a number of fixed-order perturbative computations fl(| [0], |i~2|. 

As the loop order increases, so does the computational effort. The sheer enumeration of 
various diagrams and their symmetry factors becomes non-trivial. The group theoretic and 
Lorentz structures of single graphs are involved. Finally, the scalar integrals remaining are 
hard to evaluate analytically. It is therefore clear that ideally, one would like to automatise 
the whole procedure (for a review of the current status see, e.g., |l3|| ). 

In this paper we concentrate on the first step of any multiloop computation, the enumer- 
ation of various Feynman diagrams. This step should be the easiest to automatise, since all 
one needs is a straightforward evaluation of Wick contractions. Indeed, various packages, 
such as FeynArts [14] and QGRAF [15], are available for determining n-point functions in a 
given particle physics model. 

For vacuum graphs in condensed matter systems a similar approach is possible. For the 
quartic O(N) scalar model the combinatorics is not yet too hard, but variants thereof already 
require some work. Consequently, graphical algorithms have been developed at 4-loop order 
and beyond for a number of simple models fl~6[| . 

In many cases, though, a straightforward generation of the full set of diagrams of a given 
loop order may not be the ideal way to go. In realistic theories there are very many graphs, 
and all integrals would have to be evaluated on the same footing. This is almost impossible, 
particularly if many different masses appear. 

Here we wish to present what would seem to us to be a maximally manageable setup. All 
vacuum graphs are generated, but they are cleanly separated into two groups: one, of two- 
particle-irreducible (2PI) "skeletons" with free propagators; the other, of "ring diagrams" 
with various self-energy insertions (see also (IT]]). The self-energies, in turn, are directly 
obtained from lower order skeletons. We find that this setup economises the generation of 
the various graphs quite significantly. We also point out that in low dimensions, relevant for 
statistical physics applications, the integrations remaining are qualitatively different in the 
two sets. 
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As an illustration of the setup, we enumerate the diagrams contributing to the 4-loop 
free energy of finite temperature QCD (as well as QED and the symmetric phases of the 
electroweak theory and scalar electrodynamics). We hope, though, that the setup may be 
applicable to some other cases as well. That is why we wish to separate it from the evaluation 
of the integrals arising in the finite T context [18|, specific for that physical situation. 

Our plan is the following. We summarise our basic notation in Sec. ^, reorganise the 
standard skeleton expansion in Sec. |3|, review the Schwinger-Dyson equations for n-point and 
vacuum graphs in Sec. |3], and combine them with the modified skeleton expansion to obtain 
a generating formula for skeleton diagrams in Sec. || The corresponding results are given 
for a lattice regularization of a generic model in Sec. ||. As an illustration, we show the loop 
expansion for the free energy of QCD and related models in Sec. fj]. We discuss some basic 
properties of our setup and conclude in Sec. |§[ 



2. Notation 

Let us start by introducing a concise notation. While the method is valid for any theory, we 
explicitly give all equations for a generic c^ 3 + ip 4 model. Later on we discuss more specific 
examples within this class, in particular the QCD, as well as some extensions of this class. 
The generic class also includes the electroweak sector of the Standard Model, both in its 
symmetric and its spontaneously broken phase. 
The partition function is defined as 

Z[J}= J Vipe s[v]+Jv , (1) 

where S[(p] is the action, 

S[<P] = ~\ <Pi&ij 1( Pj + ^ HjWWjVk + ^ JijkmPjWPl , (2) 

and summations over various indices, numbering (real scalar) fields and their internal and 
spacetime structures, are implied. Two comments are in order. First, we will for the moment 
not display fermions explicitly. As far as vacuum graphs are concerned, they do not introduce 
any complications apart from the usual overall minus sign for each closed loop, and can thus 
be introduced only at the end |jT9f] . Second, one should notice that the sign conventions 
in Eqs. ([[]), (|2|) are such that in the case of Euclidean actions, jijki is typically negative. 

For a theory with a broken symmetry, the inverse free propagator A -1 and the couplings 
lijk... are functions of the order parameter, but otherwise there are no essential complications. 
We return to this point in Sec. |4.l| . 

The partition function Z[J] in Eq. (jl]) is the generating functional for full Green's functions, 
r^ u11 = 5 1 } Z[J]\ J=Q . As usual we define 

W[J]=1nZ[J], (3) 

the generating functional of connected Green's functions, r™ nn = <5™ VF[</]|j =0 - Finally, one 
can define the effective action via 

S eS [4>] =W[J]- <PJ, <j> = SjW[J], (4) 
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which generates one-particle-irreducible (1PI) Green's functions, T^ PI = S e s[4>]\ ,/,=()■ Note, 
in particular, that d^Scs^] = — J. The vacuum, or free energy F (made dimensionless by a 
division with the temperature T), can be obtained from any of the generating functionals as 

F = -]nZ[0] = -W[0] = -5eff[0]. (5) 

From the basic relations <\> = 5 jW[J], ^Seff^] = — J it follows that 

t 2 J W[J]tf,S eS [<f>] = -l. (6) 

Defining, as usual, the "proper" self-energy by 

^effM^-A-i + II, (7) 

we see from Eq. @ that 5jiy[J] is the full propagator: 

S 2 jW[J] = D[(f>] = — — j- = A + AnA + AnAnA + ... . (8) 

We shall use here the following notation for free and full propagators, the proper self-energy, 
as well as general 1PI vertices: 

A = (free propagator), (9) 

D = — = + © + @ — ©— +... (full propagator), (10) 

II = @ (proper self-energy, with legs "amputated"), (11) 

A -1 = f (inverse free propagator, with legs amputated), (12) 

S^S^s = - 90 (general amputated 1PI vertex) . (13) 



n 



3. Skeleton expansion with free propagators 

We next review the skeleton expansion for the free energy F [20, 21], and modify it such that 



full propagators can be replaced with free propagators 17]. By a skeleton we mean a 2PI 



vacuum diagram: one that remains connected even if any two lines are cut. The skeleton 



expansion has been used as the starting point also in [17]. 



It can be shown [2C, 21] that the loop expansion for Eq. (|E|) can be written as 

F[D] = J2°i (Trlni^r 1 + TrlL^A) - $[D] , (14) 



3 



where i = {bosons, fermions}, Cboson = 1/2, and Cf er mion = —I- Here &[D] collects all 2PI 
vacuum diagrams. The full propagators D\ are related to their corresponding self-energies 
by D^ 1 = A -1 — II (cf. Eq. @), where A are the free propagators. Both F, U and <& can 
be regarded as functionals of the full propagators. The partition function has an extremal 
property, such that the variation of F with respect to any of the full propagators vanishes [p0| , 
[2l], |2^], giving a relation between skeletons and self-energies: 



5 D .$[D] = Ci U[D] 



(15) 



Here, we have introduced the implicit notation that whenever a term is multiplied by Cj, 
the n's and -D's following it are assumed to carry the same subscript. Pictorially, Eq. ( |i~5| ) 
corresponds to getting a self-energy by "cutting a propagator" in all possible ways in the set 
of vacuum skeletons. Hence, knowing the skeletons alone provides full information. 

In Eqs. (|I4|), (|l5|), it is the full propagators D which appear in the skeleton graphs and self- 
energies. We would instead like to obtain skeletons with free propagators. As a first step in 
this direction, we expand D in terms of the self-energy insertions II [D], D = A^ n >o(HA) n , 
to get 



F 



£ Ci Tr 



InA- 1 + (HA)' 



n>2 



*[A£(IIA)' 



(16) 



n>0 



We then have to evaluate n[D]. 

To go forward more explicitly, we restrict ourselves to 5-loop level here. Let the subscript 
n denote the loop order, and write n = Xm>i H n - It turns out that we need at most In 
a straightforward way, we obtain 



IIi 

n 2 
n 3 



HxIA^nHA], 

n^fA] + (ni rr [A + AnfA] 
nr[A] + n r 2 cd(1) [A], 

n^fA] + (n| r [A + Ani rr A]) + (ni rr [A + AnA + AnfAnfA] 



= n^ rr [A] + n r 3 cd{1) [A] + n r 3 eQ ^ [A] , 



red(2) , 



(17) 
(18) 
(19) 



where W™ are n-loop 1PI graphs, while iln ed ^ are obtained by cutting m lines in a lower 
order n^ r [A], and dressing them appropriately: 



n 



red(l) 



n^ ed(1) [A] 



n 



red(2) 



[A] = (AnrA)j5 A .nf[A], (20) 
(An 1 1 rr A) j <5 Aj n 1 2 rr [A] + (An 2 A + Anf Anf A) j( j Aj nf [A], (21) 

[A] = ~ (An i 1 rr A)j(An 1 1 rr A) k (5 A .(5 Ak n i 1 rr [A] . (22) 



For the explicit diagrammatic characteristics of IL^^, see Sec. 5.2. 
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It is easy now to unfold the loop expansion also for <&[D] = J2 n> 2^n, the last term in 
Eq. (|l6|). Up to 5-loop level, we can write 



$ 2 [L>] 



n<5 

n<5 
n<5 
n<5 



($> 2 [A + A(ni + n 2 + n 3 )A 

+A(n x + n 2 )A(iii + n 2 )A + AriiAriiAiiiA 

U 3 [A + A(n x + n 2 )A + AriiAiiiA]) 

V / n<5 

($> 4 [A + AniA]) 



n<5 



i<5 



$5 [A] , 



(23) 
(24) 
(25) 
(26) 



where the arguments are to be Taylor expanded, with first derivatives obeying (cf. the dia- 
grammatic identity Eq. (^), evaluated with free propagators) 



5 Ai * n [A]=c i lE?_ 1 [A], 



and higher ones bringing back reducible self-energies, defined in Eqs. (^0)— (22). 
Inserting these expansions into Eq. (|l~6| ) we finally get, up to 5-loop level, 



F = -^CiTrlnA" 1 + $ 2 [A 

i 

'1 



(An 



- (An 1 r + An 1 A(n 1 2 rr + -ir 2 ' 
l 

4 



red(l) 



+ $ 3 [A] + 5>iTr 

i 

+ * 4 [A]+^CiTr 

i 

+ <D 5 [A]+]T Ci Tr[i (An 1 ) 4 + (An 1 ) 2 A(nr + ^n r 2 ed{1 



i 



+ - An' rr A n' rr + n 



r red(l) 



+ Arii a 



(n|j 



red(l) 



red(2) 



or, written diagrammatically (and denoting by Fq the non-interacting result), 
— F = -F + $ 2 [A] 



+ ^ 3 [A]+^ Ci (*<o0) 

+ (^[A] + (i(0> + O + ^ O ) 



+ 



(27) 



(28) 



(29) 



5 



Here a circle with n inside denotes H 1 ", a square ifn , and a double square UT (2) . We 
will term the skeletons with free propagators, <J? n [A], irreducible. Note that the numerical 
factors in front of various types of ring diagrams do not appear to trivially follow from any 
simple symmetry argument (particularly in the case of reducible self-energy insertions), but 
are best worked out explicitly via the Taylor expansions we have described. 

Eq. (|29] ) is the starting point of our setup. It expresses the free energy in an economic way 
in terms of the irreducible skeletons $ n [A]: either as direct contributions, or as self-energy 
insertions, obtained from the same skeletons via Eqs. (|27|), (p0|)-([2^). We note that at n-loop 
level, one needs 3> n [A] but only II n _2[A], obtained from «3> n _i[A]. 



4. Schwinger-Dyson equations with full propagators 

Next, we need to generate the skeletons <3? n [A], needed in Sec. ||. To do that, we first review 
briefly the general setup of Schwinger-Dyson (SD) equations, converted to our notation. The 
SD equations will then play a central role in our main result, Eq. ([49]), which is an explicit 
formula allowing for a systematic generation of all skeletons 3> n [A] - in principle to any 
order. In this section, we follow closely the very enjoyable presentation by Cvitanovic [19]. 

4.1. General n-point functions 

The basic SD equation for the generating functional Z[J] of full Green's functions, derives 
from the trivial fact that the integral of a total derivative vanishes: 

= Jv<p5^e s ^ +J ^ = (S'[6j] + J) Z[J] . (30) 

For the generating functional of the connected Green's functions, Eq. (|3|), one gets 

= S'[W'[J]+5j} + J. (31) 

Finally, for the effective action, Eq. (||), we use from Sec. ^ that W'[J] = (ft, 5j = {8(ft/8J)8^ = 
W"[J]8$ = D[4>}6^, and J = -S'^ty] to obtain 

S' eS [<f>} = S'te + Dlfflt]. (32) 

Putting (j) — > on the right-hand-side, this gives the SD equation for the 1-point function, 
while taking derivatives with respect to <f> on both sides of Eq. (^) and putting — > only 
afterwards, generates SD equations for higher-point Green's functions, 

rf 1 = d^S'^ + DW^. (33) 

Here D[(f>] is in Eq. @, and we note that 

5^} = D[<p] (5lS cS [<p}) D^}. (34) 

A note may be in order here, concerning theories with spontaneously broken symmetries. 
In that case, (ft corresponds to the fluctuating field around some reference value v, typically 
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v = ((p). The quantity we should ultimately be computing is the free energy density as a 
function of v, i.e. the effective potential V(v) = ^/(volume). Then everything goes as before: 
we still put <j> — ► in the equations above after differentiation, while the condensate v appears 
as a parameter in the free propagators as well as in the cubic and quartic couplings in Eq. @ 
(the term J if linear in ip in Eq. (||) need not be changed p3|1). The graphs also remain the 
same: only 1PI graphs, generated by the loop expansion in Eq. (p9|), are to be included p3|| . 
Tadpole type graphs often associated with broken symmetries would only be generated if 
we want to re-expand the value of V(v) at the broken minimum in a strict loop expansion: 
writing V = En>o v n, « = En>o v n, such that Vq(v ) = 0, implies 



V(v)\ 



V'(v)=0 



+ 



+ 



vb(u ) + FiM + 



V-2 



1 {Vj) 2 

2 y " 



v 3 



v 4 



V" 
l 'o 



+ 



V=VQ 

1 (V{) 2 Vf 1 (V{) 3 V£" 
6 



- (W 

1 (Vtf + 2V{V£ 



K) 3 



+ 



1 2V{VjV{' + (V{) 2 V£ 



V=V() 



2 y " ' 2 (f ") 2 
i 3(y 1 , ) 2 (y 1 // ) 2 + 3(y{) 2 vw + {v(fv{' 

6 (F ") 3 
1 12{V() 3 V{%" + {VlfVg" 



i (v{)*(y m *- 



24 



+ ... , (35) 



where the latter terms inside the square brackets correspond to various tadpole graphs, with 
obvious notation: 1/Vq is the free propagator of the Higgs particle with a vanishing momen- 
tum, V{ iyl') is a 1-loop diagram with one leg (two legs), Vq" is a three- vertex, etc. 

Let us now illustrate the structure of Eq. ( p3| ) for the generic model in Eq. (J2]). Starting 
from Eq. @, writing down indices, and employing Eq. (|34]), we obtain for the right-hand-side 
ofEq. pi), 



1 



^k + Djkj 



+ ^Iijki{4>j4>k4>i + D jk <j)i + D M <f>j + Dijtpk + D jm D kn Di 5^ m 5^J^ o S c fi[(t)\}. (36) 

We now take further derivatives according to Eq. (|33[). Putting = after each differentia- 
tion, we thus obtain the standard equations (written in the notation of Eqs. (p|)- (|l3|) ) 



o 



(37) 
(38) 



- -->---+ --©--, (39) 

i % iii 

4- ,Q, + iO\ + K ..Q.. + .Gi. + .0, +cyciic(2 ' 3)N 



I * * 

• O (40) 

X = x • 0000 • • ' ; • •: o 

+ {2-loop terms} , (41) 

where "cyclic(ni, ri2, ...)" denotes cyclic permutations of the legs numbered. We have not 
written down the 2-loop terms in Eq. (|4l|), since they are not needed in our explicit 4-loop 
demonstration below. Likewise, all higher point 1PI functions T^ PI , n > 5, start with 1-loop 
graphs in the model of Eq. (||), and will again not contribute at this order; they will for $5, 
as well as in the model of Sec. |6[ 

Let us stress that in a local theory the manipulations needed in Eq. ([33]) can essentially be 
made using regular derivatives, and can thus easily be implemented algebraically. Introduc- 
ing furthermore h as a loop counting parameter |24[ |, allows for an iterative solution of the 
corresponding SD equations. 



4.2. Vacuum diagrams 

The SD formalism above provides equations relating n-point Green's functions. To incor- 
porate vacuum diagrams, one can use another simple trick: scaling. Noting that, e.g., Z[J] 
is a functional of all interaction parameters present in the action, Z[J, 7^, 7«&, . . .], one can 
derive hosts of relations by varying any of these parameters. 

A most useful example is to rescale the entire action as S[(p] — > ^ S[ip], and then vary h: 

-hd n lnZ[J] = Q:S[<p]) = ^ls[8j]Z[J]. (42) 

Rewriting this in the "connected" language (recall W = In Z) , 

-hd h W[J] = ±S\W'[J]+6j], (43) 

allows to finally go over to 1PI functions {d^W = d%S e R + S' eS d%<p + Jd%<j) = 5^5 e fj, W = (j) 
and 8 j = W"5<p = D[(f)}8^): 

-hd n S eS [<j>) = (~S[<p]\ = ^S[<f> + D[^}. (44) 

The free energy F = — S e s[0] can now be obtained by setting = and integrating over h. 

Noting again that after a rescaling of the integration variables, an expansion in fi is equiva- 
lent to the loop expansion one can integrate the left-hand-side of Eq. fl44|) by J n (l/h)[...], 
but on the right-hand-side one integrates over the loop number. Writing 

-S cS [0] = F = F + F mt = F + J2K nt , (45) 

n>2 
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where n counts the number of loops, it follows that 



F 



int. 



n 



n > 2 . 



Illustrating Eq. ( |46| ) for our generic theory in Eq. @, we get 



pint 



n 



(46) 



(47) 



where we again use the notation of Eqs. (p[)-(|l^). 

In principle the whole loop expansion can now be generated from Eq. (fT?|), using Eqs. fl38| ) 
The n-loop vacuum diagrams are expressed in terms of 1PI n-point functions, which in 
turn are governed by a set of SD equations. Looking closer at it though, it is somewhat of a 
mess: one has to expand full propagators in terms of free ones and the IPs, use SD equations 
to iterate loops for ITs, which brings back full propagators, etc. Fortunately, none of this is 



necessary for Eq. (29), as we now explain. 



5. Generating the irreducible skeletons $[A] 

The key observation for combining Schwinger-Dyson equations and the skeleton notation in 
a useful way is that we need to extract from Eq. ( |47| ) only a specific part, $[A]: we already 
know, by Eq. (|29|) , what all the rest combines into. But then, full propagators can be replaced 
by free propagators in all but the first term in Eq. (^7|)! Indeed, any self-energy insertion 
within one of the other graphs leads to a two-particle-reducible (2PR) diagram. For the same 
reason, the 1PI vertices in Eq. (fl7| ) can be iterated by using the SD equations of the form 



in Eqs. ([40|) , (|4T|) , but with free propagators! More precisely, it goes as follows. 

To generate the irreducible skeletons [A] from Eq. fl47| ) , it is sufficient to expand the first 
term as 



Q 



+ ( K J+ {2PR} 
Tr l+i Q + 1 QO +h (J) +\ Q) + {2PR}, (48) 



where in the second step Eq. (|3^) was used. Taking into account the minus sign in the 
relation of F and 3>[A], cf. Eq. (|29|), and writing again the loop expansion as <3? = J2 n >2 
one finally obtains a closed exact equation: 



Eq. ( |49| ) is our main result. It generates all skeletons of all orders in the theory of Eq. (J2|) , 
once Eqs. (00), (HT1) are used (with free propagators). The skeletons, in turn, generate self- 



energies via Eq. (|27| ) and the analogues of Eqs. (|20[)-(|22|). Inserted finally into Eq. (26), we 
obtain the free energy F. 



9 



5.1. Vacuum skeletons up to 5-loop level 



The procedure of working out Eq. (|49|) is simple and mechanical and can, at least up to 4-loop 
level, even be carried out by hand, as we shall demonstrate. The only complication arising 
is the identification of equivalent topologies: the same graph can be written in very many 
different ways. In order to deal with this situation, it appears easiest to assign an algebraic 
notation for the different topologies, rather than a mere graphical one. For example, one can 
count the numbers of 3-point and 4-point vertices appearing in the graph, and within those 
equivalence classes, one can use a matrix notation for how the vertices are connected. The 
significant entries of the matrix can be ordered to a single number, and by doing the same 
for all possible orderings of the vertices, a unique representative (say, the smallest of such 
numbers) can be assigned to each topology. For an explicit implementation of this kind of a 
procedure, see the 2nd paper in |l6| . 

Let us now explicitly work out the diagram classes in Eq. (|49|) up to 4-loop level. For the 
first one, inserting Eq. ( |40| ) gives either a 2-loop graph, or 3-loop graphs to be iterated further 
on, or directly 4-loop graphs: 



Here the further iterations give 



+- 



(51) 
,(52) 
(53) 



We have dropped 5-point functions each time they appear, since in the model of Eq. (pi), they 
start with a 1-loop term, so that diagrams containing them generate higher loop orders. 

The 2nd class in Eq. ( f49| ) only contributes to $2 [A], and is trivial. For the 3rd class 
in Eq. <M), 



For the 4th class, we only need the 1-loop terms in Eq. (|4l|), 

(M ~ (£>©«(© 



(54) 



(55) 



Collecting finally these different contributions together with coefficients according to Eq. ( [49|) , 

we get 



*=*0<X> 
*3=A@ + !QD+i(ID. 



(56) 
(57) 
.(58) 
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Proceeding to higher loop orders, an automatised treatment proves essential, for the rea- 
sons outlined above. Implementing our generic formulae as well as an ordering algorithm 
separating topologies in FORM |25[], we obtain in a straightforward way the complete set of 
5-loop skeletons, 



•■= ^(H>i(H>KH>i(D + '© + '© + '(B 



© + ™SB + A(S)^(M) + iO + ™ 



(59) 



Note once more that these skeletons are all that is needed for generating the loop expansion 
for the full free energy, as discussed above. 

5.2. Self-energies up to 2-loop level 

Now that we have $ n [A] in Eqs. (|56|)-(|59|), irreducible as well as reducible self-energies can 
easily be obtained with Eqs. (|27|), (|20|)-(p2|), etc. For bosonic particles, for instance (c; = ^ ), 
we get 



n 



T-rirr 
LL 2 



red(l) 



-©- 



-©- 



o ->o 



1 - 




- +T 



(60) 
(61) 
(62) 



etc. Note that the outcome of the derivative in Eq. fl27|) must be symmetric in all (bosonic) 
indices. The 3- and 4-loop self-energies could be derived from $4 and $5, respectively, but 
we choose not to give them here, since they are not needed for the set of 4-loop vacuum 
diagrams that we will display explicitly in Sec. |7| 

With Eqs. (|60|)-(|6^), the ring diagrams in Eq. (^) are readily written down. 



6. Generic model on the lattice 

So far we have considered the generic model in Eq. (||). However, in a lattice regularization 
of gauge theories, higher vertices appear as well, without spoiling renormalizability. At the 
generic level, it is straightforward to add such couplings to the theory in Eq. (|2|). We can 
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include, e.g., terms up to ~ (l/&l) r yijklmnop t Pi ( Pj ( Pk'Pl i Pm ( Pn ( Pa l Pp] as would arise in lattice 
perturbation theory for SU(iV) gauge theories, if one keeps terms contributing to 4-loop 
vacuum graphs. Such computations would be needed when one converts results of three- 



dimensional numerical Monte Carlo studies from lattice to continuum regularization [1C]. 

In this case, everything goes as before, except for the appearance of extra vertices in the 
SD equations, as well as in Eq. (|4^). We shall here simply spell out the final results, without 
rewriting explicitly the modified SD equations. We obtain the following additional skeletons, 



$ 3 
$4 



lat 



lat 




+±©0*00^000* 




+ 



hi 





^384 




as well as the additional irreducible self-energy, 



nf 



lat 



--©-- 

where we again assumed q 



lat 



(63) 



(64) 



(65) 



7. Applications: QCD, QED, SQED, electroweak theory 

As an application of the generic formulae derived above, we consider in this section SU(iV) 
gauge theory with fermions and a scalar field. This class includes QCD and QED (where 
graphs containing scalar propagators and, in the latter case, gauge field self- interactions, are 
to be dropped out), as well as the electroweak theory and scalar electrodynamics (SQED). 
For brevity, we display here only the vertices appearing in the symmetric phases of the latter 
theories. We mostly use the language of QCD, referring to the gauge fields as gluons, etc. 

The Lagrangian is specified by giving Feynman rules for the free propagators and free 
vertices, 

— , - , — • A • A. • A • X • X • X • m > 

where gluons/scalars are denoted by wavy /straight lines. Both quarks and ghosts are denoted 
here by dotted lines; the Feynman rules for them are different, but the symmetry factors 
agree — the only exception being diagrams with more than one closed fermion loop, in which 
case both ghosts and quarks can appear in the same diagram simultaneously, reducing the 
symmetry by an obvious factor. 

We do not here write down counterterms explicitly. Coupling constant counterterms can 
be viewed as a part of the cubic and quartic couplings, while wave function and mass coun- 
terterms can be treated as a part of the irreducible self-energies IIJ£ r , making their appearance 
only in ring diagrams according to Eq. (|29|), 
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Let us first note that once we write down the summation over the field content explicitly 
in Eq. @, the "natural" symmetry factors in front of the vertices change. For instance, 
writing the 4-point vertex in the case of two sets of fields, — > {Ai} + {B a }, and using 
the symmetry of jijki, one gets 

111 1 

-^flijkmVjVWi = -^HjkiAiAjAkAi + —ji jka AiAjA k B a + —-^"f ijo ,/3AiAjB a Bp + ... . (67) 

Similarly, writing the 3-point vertex for three different fields, {(pi} — > {Ai} + {B a } + {Cm}, 
one finds 



1 



1 



1 



Qlijk t Pi t Pj t Pk = QlijkAiAjAk + —Jij a AiAiB a + ^iaMAiBaPu + ••• • 



(68) 



With these conventions, each tree-level vertex in the graphical notation corresponds just to 
JijkhJijka, etc, without any symmetry factors there: all of them are shown explicitly. 

The only thing remaining is to write the summation over particle species explicitly also in 
the propagators of Eqs. (|56j) (|58|), 



+ 



+ 



+ 



(69) 



Only the vertices allowed by the Feynman rules are kept after this substitution. This generates 
all the graphs, with the correct symmetry factors. 

7.1. Vacuum skeletons up to 4- loop level 



The procedure outlined above can easily be carried out explicitly, and up to 4-loop level even 
by hand. The main complication is again the identification of various equivalent topologies, 
and for this a suitable algebraic notation may be more useful than a graphical one. As a 
result, for the field content in Eq. (|6q), we finally obtain 



<J> 2 

$ 3 



*GO + ^©~^U + ^OQ + *0 + ^00 



21 



+ 7 



+- 



+- 



+ 



+ 



(70) 



(71) 



4—79 



-1 T 



.1/ U } . iy 
3 \ i 1 !T|! ' 



2 V ? ~ 8 



(EH 



+ 



16 



! + 



18 



(H)+*(HH(H) 



13 



+^ (TV) +^ (TT) +\(Th +i 



+ 
+ 



^4 



+: 



+1 



+ - 



+1 



+ 7 



t+7 



+ 



+ 



16 



+ 



+ 



16 



+ 



18 



(72) 



7.2. Self-energies up to 2-loop level 



Using Eqs. fl27| ), (|20f) , the skeletons above immediately produce the self-energies of the model 
in Eq. (|66|), We obtain 



(73) 
(74) 
(75) 



-©- 



+ i -^^y* + i +i -^^y +\ "C^" 

1 CD - <i> ' 
+ i <x>+i^0-+i-C> + ^O- + ^O- ' (78) 



(76) 
(77) 





(79) 
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7.3. Ring diagrams up to 4-loop level 

To be exhaustive up to 4-loop level, let us finally give the set of ring diagrams for the model 
of Eq. (|66|). While there are no ring diagrams up to 2-loop level, from Eq. (|2^ ) we get 




Note the extremely economic structure of the skeleton expansion of Eq. (|29|): the few ring 
diagrams above summarize 22 (276) 3-loop (4-loop) diagrams. 

8. Discussion 

In this paper we have described a simple practical procedure for systematically generating 
all vacuum diagrams of a given loop order in a generic field theory. 

We have shown that the sum of vacuum diagrams can be written in the form of a modified 
skeleton expansion, Eq. (p9|). It contains 2-particle-irreducible "skeletons" with free propa- 
gators, as well as various self-energy insertions inside "ring diagrams". The self-energies are, 
in turn, determined by the skeletons. Thus, all one really needs is the skeletons. 

The 2-particle-irreducible skeletons of a given order are, then, generated by Eq. ([49|). It 
contains a number of full 3-point and 4-point vertices, which can in turn be expanded using 
specific "irreducible" Schwinger-Dyson equations (Eqs. (|4C|), (|4l"|), etc), where full propagators 
have been replaced with free propagators. In this way, all vacuum graphs are generated 
simultaneously, with the correct symmetry factors. Finally, the precise particle content of 
the theory one is interested in can be specified as discussed in Sec. |7[ Our method is also 
directly applicable to theories with spontaneous symmetry breaking, as only free propagators 
and vertices are modified; tadpole graphs are generated by Eq. (|35|). 

This iterative procedure is very straightforward and can be automatised, but up to 4-loop 
level the computations are easily carried out even by hand, as we have demonstrated. Thus, 
we believe that our setup economises the generation of the set of high order vacuum diagrams, 
compared with techniques where all the types of graphs have to be dealt with on the same 
footing, without a separation into skeletons with free propagators, and ring diagrams. 

Up to this point, we have not discussed at all the integrations remaining to be carried 
out, after the diagrams have been generated. Let us end by pointing out that our setup is 
beneficial as far as their structure is considered, as well, in dimensions lower than four |Tfl| . 
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The point is that low-dimensional field theories of the type in Eq. (|2|) are super-renorma- 
lizable. In fact, for d = 2, 3, only the 2-point function suffers from ultraviolet divergences, as 
can be seen by simple power counting. Therefore the skeleton graphs which by definition do 
not have any genuine 2-point functions inside them, do not contain any ultraviolet divergences 
in subdiagrams. The ring diagrams, on the other hand, do have ultraviolet divergences in 
subdiagrams. Note, in particular, that since n^ r ,iln ed ^ m ^ come with different symmetry 
factors in Eq. (|29|), the counterterms in n^ r which make the whole n n finite, do not in 
general immediately cancel all the ultraviolet subdivergences of the ring diagrams. 

Consequently, various ring diagram classes can contribute to the overall divergences of 
the vacuum graphs with potentially infrared sensitive coefficients, coming from the other 
parts of the final integration, while skeleton diagrams can not. Fortunately, the ring diagram 
integrations are simpler than those in the skeleton graphs, and this problem can thus be dealt 
with in a tractable setting . 
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